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Abstract 

This theoretical paper presents numerical calculations for photoassociation of ultracold cesium 
atoms with a chirped laser pulse and detailed analysis of the results. In contrast with earlier 
work, the initial state is represented by a stationary continuum wavefunction. In the chosen 
example, it is shown that an important population transfer is achieved to ~ 15 vibrational levels 
in the vicinity of the v=98 bound level in the external well of the 0~(6s + 6^3/2) potential. Such 
levels lie in the energy range swept by the instantaneous frequency of the pulse, thus defining a 
"photoassociation window". Levels outside this window may be significantly excited during the 
pulse, but no population remains there after the pulse. Finally, the population transfer to the 
last vibrational levels of the ground a 3 S^(6s + 6s) is significant, making stable molecules. The 
results are interpreted in the framework of a two state model as an adiabatic inversion mechanism, 
efficient only within the photoassociation window. The large value found for the photoassociation 
rate suggests promising applications. The present chirp has been designed in view of creating a 
vibrational wavepacket in the excited state which is focussing at the barrier of the double well 
potential. 

PACS numbers: 
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I. INTRODUCTION 



After ultracold atoms, ultracold molecules are presently a subject of constant interest, 
stimulated further by the recent observation of molecular condensates . This is wh 

the various routes leading to the formation of ultracold molecules are actively explored 
Starting from molecular beams, two non-optical techniques, buffer gas cooling of molecules 



dry 

J 



now reach temperatures well below IK. 



and Stark deceleration of polar molecules [t], 0] 
Another route relies on optical techniques: laser fields are used to cool down alkali atoms, and 
to create molecules via the photoassociation reaction 9]. Short-lived molecules formed in an 
excited electronic state must then be stabilized, by spontaneous emission or other radiative 

nnnnnn 

coupling, into bound vibrational levels of the ground electronic state [1C| [11|, 112J, [13|, [15|| . 
The translational temperatures are much lower (T < 20/iK) and higher densities are to be 
expected than with non-optical techniques. An important drawback, however, is that the 
stable molecules are produced in a superposition of vibrational levels, among which some are 
very excited. Bringing such molecules to the v =0 level of the ground electronic state, thus 
reaching ultralow vibrational temperatures, is therefore an important issue. More generally, 
finding new photoassociation and stabilization schemes is an interesting research subject. 

Up to now, photoassociation experiments have mostly been using continuous lasers: a 
more active role of the radiation would be to optimize the reaction by shaping laser pulses. 
The aim of the present paper and of the following ones is to explore the feasibility and ad- 
vantages of experiments using chirped pulses to form ultracold molecules, therefore bridging 
the gap between the subjects of ultracold matter and coherent control. 



The field of coherent control has developed both theoretically 



experimentally 
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with quite different applications, such as photodissociation of 
molecules. It has indeed been demonstrated that shaping laser pulses could enhance consid- 
erably the yield of a photodissociation reaction, so that a similar effect is expected for the 
reverse reaction. 



Very few papers consider the photoassociation reaction wit 



mal regime, ultrafast photoassociation has been studied |24 , 



i pulsed lasers: in the ther- 
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271 ] , showing the validity of 



the impulsive approximation. In the ultracold regime, several theoretical and experimental 



papers have proposed a time-dependent study of the photoassociation reaction 
eventually considering also the formation of long-lived ultracold molecules |3] 
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Most the- 
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oretical calculations are using a wavepacket representation for the initial state of the two 
colliding atoms. In this framework, considering a wavepacket localized at the outer classical 
turning point of the photoassociated level, our group has discussed the characteristic times 
for photoassociation of cold cesium atoms [31] and the separation of time scales. 

Photoassociation of cold atoms with a chirped-pulsed-laser field was first theoretically 
explored by Vala et al j^ : appreciating the slow nuclear relative motion in ultracold colli- 
sions, a chirped pulse was optimized to achieve total transfer of popu lation under adiabatic 



following conditions proposed by Cao, Bardeen and Wilson 



34j . l35j | . In these calculations 



the initial state was represented by a gaussian wavepacket and a total population trans- 



fer was indeed obtained in the region o 
drawback of both this study and Ref. 



distances where this wavepacket is localized. The 
was that they overlooked the precise nature of 



the initial continuum state. Since the collision energy is so close to threshold, a plane wave 
description or a superposition of Gaussian wavepackets do not address the actual shape of 
the wavefunction. A correct representation of the initial continuum state of two ultracold 
colliding atoms should be a stationary continuum state: indeed, experiments with continu- 
ous photoassociation lasers identify several nodes of the stationary wavefunction as minima 
in the photoassociation signal 3^ . 

The aim of the present paper is to go further in that direction, considering now pho- 
toassociation with chirped laser pulses when the initial state is a delocalized continuum 
wavefunction representing the relative motion of two atoms at very low collision energies 
kT. Such a wavefunction displays very slow oscillations at large distances, and numerical 
calculations must then use large spatial grids. On the other hand, for ultracold processes, 
where the initial continuum state is very close to threshold, even a very weak light field 
couples the last bound vibrational levels of the ground electronic surface, which must also 
be correctly represented on the grid. Then, the numerical calculations need a mapping 
procedure in order to reduce the number of grid points 

It should be appreciated that the vibrational levels of the photoassociated molecule, being 
close to the dissociation limit, are physically very different from lower lying vibrational 
states. Their vibrational periods are orders of magnitude longer and most of the amplitude 
is concentrated in the outer turning point. As a result, for even weak light fields, the Rabi 
period associated with population transfer to the excited state can be shorter than the typical 
vibrational periods. For short enough pulses, the relative motion of the two atoms during 
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the pulse duration can be neglected, thus allowing simple interpretations in the framework 
of the impulsive limit for light-matter interaction j^. This implies that, during the pulse, 
the wavefunction at each nuclear separation can be decoupled from other positions, so that 
each distance can be considered as a radiation coupled two-level-system, where transfer of 
population is analyzed. Another limit now emerging is the limit of adiabatic transfer, where 
the intensity of the light is sufficient so that the associated Rabi period is smaller than all 
other relevant timescales. 

The paper shows how a two-state model can be developed to provide analysis of the 
results within an adiabatic population inversion mechanism, effective only within a range 
of energies defining a photoassociation window. Also discussed is the possibility of avoiding 
population transfer outside this window. An example of application is the focussing process, 
where photoassociation is creating in the excited electronic state a vibrational wavepacket 
optimized for transfer to a vibrational level of the ground state via a second pulse. 

The paper is organized as follows: in Section |TI] we recall the theoretical model for pho- 
toassociation with chirped laser pulses. In Section IIHI we discuss the various timescales of 
the problem. In Section IIVI we give an example of numerical calculations for Cs2 0~ pho- 
toassociation, with a chosen chirped pulse, where a few neighbouring vibrational levels stay 
populated after the pulse. For interpretation, section IS develops, in the framework of the 
impulsive limit, a two-state adiabatic model for population inversion within a photoasso- 
ciation window. In Section IVI[ the validity of this adiabatic approximation is discussed. 
Section IVHI is the discussion of the numerical results in the frame of the simple model for 
adiabatic population transfer. Finally, in Section IVlIII we show how the chirp parameter can 
be chosen in view of focussing a vibrational wavepacket of the photoassociated molecules, 
such as control of the formation of stable ultracold molecules can be performed. Section IIXI 
is the conclusion. 

In a forthcoming paper, hereafter referred to as paper II, we shall give more examples 
for numerical calculations and discuss optimization of the photoassociation yield and of the 
cold molecule formation rates. 
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II. THEORETICAL MODEL 



A. The photoassociation reaction 

The photoassociation reaction starts from two cold atoms, at a temperature T, colliding 
in the ground state potential V groun d, and which absorb a photon red-detuned compared to 
a resonance line to yield a molecule in a vibrational level of the excited electronic potential 
V eX c- We shall consider the example of two cold cesium atoms colliding in the lower triplet 
ct 3 E+(6s + 6s) potential, and forming a molecule in a vibrational level of the 0~(6s + 6p 3 / 2 ) 
potential: 

Cs(6s 2 S 1/2 ) + Cs(6s 2 S 1/2 ) + hu L - Cs 2 (0- (6s 2 S 1/2 + 6p 2 P 3/2 ; v, J)). (1) 

The reaction ((TJ) is the usual representation of the photoassociation, using a continuous laser 
red-detuned by 5£* relative to the D 2 resonance line: 

hu L = hu at - sf, (2) 

where hu at is the energy of the atomic transition 6s — > 6P3/2, which excites a single 
vibrational level of the 0~ electronic state. In the present work, we shall treat excitation 
with a chirped laser pulse characterized by time-dependent frequency uj(t) or detuning 
5£*(t), and which populates several vibrational levels. 



As indicated in FigJTJ the 0~(6s + 6p 2 P 3 / 2 ) potential curve displays a double well 
structure with a hump located around 15 do- The levels populated by the photoassociation 
reaction correspond to vibrational motion in the outer well: Vtot is the vibrational quantum 
number in the full potential, while v is the numbering of the levels in the external well. 
In the experimental photoassociation spectrum, vibrational levels from v=0 to v =132 are 
identified For two levels, a tunneling effect is present which has been analyzed from 
the photoassociation spectrum 

The theoretical model that we are presenting is of course not connected to the particular 
shape of the excited potential V exc . 
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B. The two-channel coupled equations 



The dynamics in the ground state and the excited state is described by the time- dependent 
Schrodinger equation 

H*(#, t) = [H mol + W(t)MR, t) = ih^*(R, t). (3) 

where *&(R, t) is a two-component wavefunction, with ^ gr0 und(R, t) and ty exc (R, t) describing 
the radial relative motion of the nuclei in the ground and excited potential respectively. 
The molecular Hamiltonian H mo i = T + V e i is the sum of the kinetic energy operator T 
and electronic potential energy operator V el . The coupling term is written in the dipole 
approximation: 

W(t) = -£>(R) ■ ei S(t) (4) 

involving the dipole moment operator D(i?) and the electric field defined by a polarization 
vector ei (assumed to be constant) and by an amplitude S(t). We define D ge (R) from 
the matrix elements of the dipole moment operator components between the ground and 
the excited molecular electronic states. Since the photoassociation reaction occurs at large 

— * _ 

distances, we shall neglect the R— dependence, using the asymptotic value D ge (R)-el ~ Z?|£ 
deduced from standard long-range calculations 41 1 . For this reason, the formulas below are 
written for a R-independent coupling. 

1. Excitation by a continuous laser 

For a continuous laser with constant frequency ujl/2it, the electric field is 

S(t) = So cos(u it) (5) 
where So is the amplitude of the field. 

The explicit temporal dependence of the Hamiltonian H is eliminated in the framework 
of the rotating wave approximation, defining new radial wavefunctions in both excited and 
ground state channels through 

*^( J R,t) = exp(i^)*« ec (i2,t) 
^ g L (R,t) = eM-i^r)y 9 round{R,t), (6) 



and neglecting the high frequency component in the coupling term. This allows to write the 
radial coupled equations as 



in— 



( 



T + V^R) W L \ I ^ e L {R,t) 
dt \*p(R,t) ) \W L T + V' g (R) ) \^(R,t) 

In Eqs. the two non diagonal terms are identical and time-independent, 



(7) 



where I is the constant laser intensity, c the velocity of the light and eo the vacuum permi- 
tivity. 

The diagonal terms in Eqs. ((7j) involve the dressed potentials 

V e (R) = V exc (R) - hu L /2 = V{R) + A L (R), V' g {R) = V ground {R) + hu L /2 = V(R) - A L (R). 

(9) 

In Eq. (jHJ), we have introduced the difference between the two dressed potentials: 

2A L (R) = V' e {R) - V' g {R) = V exc {R) - V ground (R) - hu L , (10) 
where 2Ai{R) — > 5^ for large R, and the mean potential 

\T(J^ — Vexc(R) + Vground(R) 

When the dressed-potentials are crossing at a large distance Rl (V^(Rl) = V 9 {Rl)) 

A L (R L ) = 0, (12) 
rotoassociation into a vibrational level with binding energy 



the resonance condition for 
may be estimated through 



4] 



V e (R L ) - V e (oo) « -E%, (13) 



where we have neglected the kinetic energy of the initial continuum state. Therefore, the 
outer classical turning point of the photoassociated level should be located close to Rl. 

2. Chirped pulse and transform-limited pulse 

In order to optimize the formation rate of molecules, we shall consider gaussian "chirped" 
pulses (34), where the laser field S(t) has a quadratic time- dependent phase (p(t), hence a 



7 



time-dependent frequency u(t) together with a time-dependent amplitude So fit). We have 
to note that the discussion reported here is general and can be adapted for other types of 
temporal envelopes than the gaussian one employed here, as for more general chirp rates 
that could be defined, beyond the linear approximation (Eq.(|18|)) used in the present paper. 

• the field is 

£(t) = E f(t)cos[uj L t + <p(t)}, (14) 
where the amplitude involves a gaussian envelope fit) centered at t = tp, 

/(t) = ,/^exp[-21n2(^) 2 ], (15) 

V TC T C 
TL = 1 = < T C (16) 

^l + ( X r2/41n2)2 

with the temporal width Tc defined as the full width at half maximum (FWHM) of 
the temporal intensity profile £lf{t) 2 . 

• the phase ip{t) in Eq. (fPfj) has a second derivative equal to the linear chirp rate 



y = 

A dt 2 ' 



1 1t~ 

<P( 1 ) = 7iX(t ~ tp) 2 ~ o arccos u L t P . (17) 

2 2 t c 



• the frequency of the field, related to the derivative of the phase, varies linearly around 
the carrier frequency ul- 

u{t)=u L + ^ = u L + x{t-t P ), (18) 

Each pulse is thus defined in the temporal domain by five parameters which are the 
carrier frequency Ul, the chirp rate x, the temporal center tp, the width tq , and £ . In 
the present study of photoassociation, the carrier frequency u>l is perfectly defined from an 
atomic resonance frequency in terms of the detuning 5£* (see Eq. ^). 

Concerning the two last parameters, we must note that for \ = 0? t l — ^c, defining 
a transform limited pulse with FWHM of |/(t)| 2 equal to tl, and maximum amplitude 
of the field S - Therefore the maximum amplitude of the chirped pulse, Em = £oy~> is 
always smaller than So'- the chirp increases the time-width of the pulse while decreasing its 
maximum amplitude. The coupling in Eqs.(J7|) now becomes time-dependent and reads 



W L f(t) < W max = W L Jp- = -MU—\D%\ (19) 

V T c V T c 2 V ce 



where we have defined a peak intensity II = ^f-^o- 

In the following, we shall assume that the relative phase of the ground and excited electronic 
wavefunctions is such that Wl > 0. 

• In contrast, the chirp does not modify the width of the pulse in the frequency domain, 
which stays proportional to 1/tl- Indeed, £(o>), which is the Fourier transform of S(t), 
displays a gaussian profile with FWHM \/25uj= 4v / 21n2/rx, and a phase $(u;) which 
is a quadratic function of the frequency, leading to the definition of a linear chirp rate 
in the frequency domain 

J2a 2 2 

$" = = y r ^ , (20) 
rfcu 2 X (41n2) 2 1 J 



Consequently one has 



($")2 



T c = T 2 [l + (41n2)^]. (21) 



/ |£(t)| 2 dt = 0.98£ pu , se (23) 

J —Tn 



Finally, the chirp does not change the energy E pu i se carried by the field, which is 
proportional to the square of the amplitude So and to the temporal width tl of the 
transform limited pulse 

This allows us to define a time window, since for a gaussian pulse 98 % of the energy 
is carried during the time interval [—Tc, +Tc)' 

ce f +T c 

Since from Eq.(fTK|) f(tp±rc) = \\J^) we ma y a l so define a lower limit to the coupling 
term during the temporal window and write 

V(t - t P ) E [ Tc, +t c ] ^ 7^ L ,/^ < W r /(t) < ^l./^ = WW (24) 

4 V r c V r c 

The definition of such a time window will reveal very useful for the analysis of the 
dynamics. 

We have represented in Fig. |2k the typical variation of the amplitude of a transform limited 
pulse and of the pulse obtained through linear chirp, illustrating both the stretching of the 
temporal width from tl to tc, the reduction of the peak intensity, and the conservation of 
the energy carried by the field. The energy range swept during the time-variation of the 
frequency (see Eq.(fT%|)) is illustrated in Fig. |2b- 



3. Coupled radial equations for excitation with a chirped pulse 

When considering excitation via a chirped pulse, we may apply the rotating wave approx- 
imation as previously at the carrier frequency uj l . Once the terms oscillating as ex$[±2iuj L i\ 
are eliminated, the resulting coupling term is time-dependent with a phase variation, so that 
the coupled equations now are 



(25) 



T + V(R) + A L (R) W L f(t)exp[-Mt)] ^ L (R,t) 
W L f(t)exp[Mt)} f + V(R)-A L (R) ) \*p(R,t) 

However, the existence of a complex coupling term makes the discussion intricated: it is 
therefore convenient to modify further the rotating frame by defining new wavefunctions for 
the excited and the ground state, 



VZ(R,t) = exp(^/2) x V^(R,t) = exp( ^ L * + V e*e(fl,0, 



W!(R,t) = exp(-i<p/2) x ^ q L (R,t) = exp(-i 



00 L t + If 



)^ ground^R-i t)i 



(26) 



thus performing the rotating wave approximation with the instantaneous frequency. The 
coupled equations now read 



*" e (R,t) 



T + V{R) + A L (R)-*% W L f(t) 



V 



W L f(t) $ + V(R)-A L (R) + %%) \ *?(R,t) J 

In the following, we shall call W(t) the real time-dependent quantity 

W(t) = W L f(t), 

and 2A(R,t) the time-dependent energy difference such that 

2A(R, t) = 2A L (R) -h^ = 2A L (R) - h X (t - t P ) 



(27) 



(28) 



(29) 
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4- Resonance window and spectral width 



It is important to note that, in contrast with cw excitation, the relative position of the 
two potential curves at a given distance R now varies linearly as a function of t. Hence, at 
some distances the upper curve may become the lower or vice versa, while the position of 
their crossing point Rc(t), such that 

A(R c (t),t)=0; R c (t P )=R L (30) 

varies with time around the distance Rl defined in Eq. (I12j) . 

We shall consider situations where during the time window [tp — Tc,tp + tc], the crossing 
distance Rc(t) is spanning a range of distances 

R-min — R — Rmax 

X 

Rmin = Rc(tp — | — rTc) 2A L (R min ) = — h\X\Tc, 

\x\ 

X 

Rmax = Rcitp + J—Jc) => 2A L (R max ) = h\x\r c - (31) 
I X I 

The definition of such a resonance window requires two conditions on the width of the pulse 
and its frequency, to ensure that the curves keep crossing: 

h\ X \rc < U™ n (32) 
Mx\rc < 8? (33) 

where —U™ tn is the minimum value of the difference potential 2Ap{R) at R < Rl , while 
5ff the asymptotic splitting. The resonance window is defined independently of the pulse 
intensity. However, it depends upon our choice for the temporal window, and this point will 
be further discussed below. 

Writing the resonance condition from Eq. (J3U|) . which depends upon the central frequency 
of the pulse, would be misleading: due to the spectral width, a range of levels in the 
neighbourhood of the resonant level can be excited. We shall discuss in more details these 
two aspects of the excitation with a particular example. 
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III. TIME SCALES 



A. Time scales related to the radiation 

The two radiative phenomena in our problem are the spontaneous emission time T spont 
and the Rabi period associated to the coupling with the laser. We shall study phenomena 
at a time scale short compared to T spont . For a cw laser, the period T^ abi of the "Rabi 
oscillations" associated with the coupling Wl between the two states resonantly coupled is 

Tkabi — yy^- (34) 

Once a pulsed laser is introduced, due to the variation of the peak intensity, this time 
constant is modified. We define two similar time constants, at the maximum intensity 
(t = tp), and at the edges of the time window (t = tp ±tc), characterizing the dynamics in 
the vicinity of the instantaneous crossing points Rl, R m i n , Rmax' 

lRabAtp) ~ W max ~ir L ia * 
Tlu (t P ±r c ) = §^= AT c Rabl (t P ) (35) 

*V m.n.T. 



As a generalization, following Ref . 42] , we shall also use a local time-dependent Rabi period 
defined by 

TRab(R,t)= = 36 

^fW 2 (t) + A*(R,t) 

Obviously the pulse duration Tc is also an important characteristic of the pulse. For the 
discussion of the adiabaticity condition, and for the optimization of the pulse, it is interesting 
to evaluate the time constant associated with the energy range swept during the time window 

T c hirp TTi i (37) 

2|x|t c 

This quantity should be compared to the width of the energy distribution |£(c<j)| 2 , which 
is 5uj = 41n2/r^, and does not depend on the chirp rate. A time scale caracteristic of the 
spectral width is then introduced as 

2tt 2tt 

T SV ect = t- = -; — n = 2.26r L . 38 
p 5u 4 n2 v ; 
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The ratio p between those two time constants is bounded, since 



p = IfE^L = 2 l*l r c = 2x (39) 

Tchirp 5uJ y/1 + X 2 

< p < 2. (41) 

Therefore the frequency band swept by the central frequency due to the chirp is never 
larger than twice the spectral width. This is an important parameter in the description of 
the chirp process. 



B. Time constants associated with the dynamics. 

When discussing the vibrational motion of a vibrational level v, with binding energy E v 
in the excited state, we shall consider the classical vibrational period estimated from 

T mb (v) ~ 2^|| « ^ . (42) 

For the discussion on optimization of the chirped pulse, it is also relevant to consider the 
revival period, defined |43J as 

m , s 47r/i 

TreM « — — ■ (43) 

For the ground state, we shall consider the motion at a low collision energy E, described by 
a scattering length L corresponding to a time-delay 

J S (E) HL /AA . 

T = h -glT ~ Ik ' (44) 

where we have introduced the elastic phaseshift 5(E) = Arctg(-kL), the reduced mass p 
and the wavenumber k. 



IV. AN EXAMPLE OF NUMERICAL CALCULATIONS FOR PHOTOASSOCIA- 
TION WITH A CHIRPED PULSE. 

Since the pulse depends upon many parameters, the way of optimizing them is not 
straightforward. We give results of numerical calculations in the example already presented 
in Section III Al of photoassociation into several bound levels in the external well of the Cs2 
g potential, populated by a chirped pulse. 
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A. The physical problem and its time scales 

The photoassociation reaction uses a chirped laser pulse linearly polarized, already illus- 
trated in Fig|2l for which the parameters are reported in Table [Q 



8f 


II 


W L 




TL 








2.656 cm -1 
1.21 xlO" 5 au 


120 kW cm" 2 


0.7396 cm" 1 
3.37 xlO- 6 au 


0.98 cm" 1 


15 ps 


34.8 ps 


-0.025 cm" 1 ps" 1 
X=-0.28 x 10~ n au 
X=-4.79x 10~ 3 ps~ 2 


-170 ps 2 



TABLE I: Parameters for the pulse excitation considered in the present paper (see section Til B 21 
and Fig. EJ): detuning <5£* , intensity II, coupling Wl, energy associated to the spectral width hSu, 
temporal widths tl and tc, linear chirp parameters x an d 

The detuning b a L l of the carrier frequency corresponds to resonant excitation at 
t — tp of the level t>=98 of the Cs2 0~(6s + 6^3/2) potential, and to a crossing point 
Rl = 93.7 a . Calculations have also been performed by changing the sign of the 
chirp. For a continuous laser excitation with tc polarization between the electronic states 
3 E+(6s,6s) and 0~(6s + 6^3/2), and neglecting the i?-variation of the dipole coupling, the 



intensity is related to the coupling Wl by Eq. (jSJ), giving Wl(o,u) = 9.74 x 10 9 ' J iw/em 2 ill 



The time scales associated with the radiation are reported in Table |H] They appear to 
be of the same order of magnitude: in particular, T c i lirp and T spect are similar, since for the 
present pulse, the resonance width 2|x|r c is nearly twice the spectral width 5uj, leading to a 
value p = 1.8 close to the maximum value which is 2. 

The lifetime of the atomic level Cs(6p 2 P3/2) being 30.462 ns, the radiative lifetime of 
the photoassociated levels in the external well of the 0~(6s + 6^3/2) curve is ~ 30 ns. In 
the present calculations where spontaneous emission is not introduced, we shall study the 
evolution of the system at much shorter time scales. Through incertainty relations, it 
means that the definition of the energy is larger than the natural linewidth of the cesium 
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Tspont 


rpL 

1 Rabi 




T%abi(tP±Tc) 


TRab(00,tp) 


TRab{R = I5a ,tp) 


T c hirp 




P 



30 ns 



22.5 ps 



34.3 ps 



137 ps 



12 ps 



0.18 ps 



18.9 ps 



33.9 ps 



1. 



TABLE II: Time scales associated with the radiative coupling: spontaneous emission time T spont i 
Rabi period Tj^ abi , chirped Rabi periods T^ abi (tp) and T^ abi (tp ± Tc) defined in Eqs.( I55)l . 
Besides, we have reported the time-dependent local Rabi period Tp a ^ at the maximum tp of the 
pulse for R = 15 ao and R — > oo. Finally, the time constant characteristic of the spectral width of 
the pulse T spect is reported together with the chirped characteristic time T c hi rp as defined in Eqs. 
(|37|38|h p is the ratio between the resonance width and the spectral width defined in Eq.([39|) 

atom resonance line AE= 1.7 10~ 4 cm~ 1 = 5.3 MHz. The hyperfine structure splitting of 
9193 MHz will not be considered in the present paper, but it has a significant value at the 
time scale of the problem. 



At resonance (A = 0) and for a cw laser, the intensity of the pulse would correspond to a 
Rabi period of 22.5 ps. However, due to the streching factor tq/tl reducing the maximum 
intensity, it is increased to 34.3 ps at t = tp. In the wings of the pulse, this time constant is 
further increased, reaching the value 137 ps at the times tp ± tq- Such characteristic times 
should be compared to the times scales associated to the vibrational motion, and to the 
collision time. 

For the excited potential V exc (R), the outer well in the 0~(6s + 6p 3 / 2 ) curve was fitted to 
photoassociation spectra by Amiot et al ^J. The latter was obtained Q| by matching to 
ab initio calculations at short and intermediate range. The lower triplet state potential 
3 S+(6s,6s) has been chosen in order to ensure a correct value of the scattering length L, 
lere taken as L ~ 525 a : the asymptotic behaviour being Cq/R 6 , where C$= 6828 au 
44| . the short range part extracted from had to be slightly modified j3] for that purpose. 

Under the conditions illustrated by the Table HI the time-dependent energy difference 
A(R,t) between the two dressed potentials, denned in Eq. (|2~9"j). is drawn in FigJH 
The long range splitting between the two potential curves is 2Ai(oo) = <5£*=2.6 cm -1 , 
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resulting at t = tp into a local Rabi period of 12 ps in the asymptotic region, as reported 
in Table [H] The energy difference reaches large values at short distances, yielding a 
very small local Rabi period of 0.18 ps at R— 15 ao The instantaneous crossing point 
(defined by A(R,t) = 0) varies with time, and we have indicated in the Figure the 
values R min =85, Rl=93.7 and R max =107A a defined above in Eqs. (J3T|) and discussed 
in Section III B 41 Such distances are close to the outer turning points of the vibrational 
levels v=92, 98 and 106 respectively, for which Table IIHI displays parameters such as 
binding energies and vibrational periods. The pulse with negative chirp value described 
in Table IJ has an intantaneous frequency resonant with the v=106 level at t = tp — tc, 
with v=98 at t = tp, and v=92 at t = tp + Tc- Alternatively, when the chirp becomes 
positive, the resonant condition is verified first by v =92 and finally by the upper level v =106. 



Vtot 


V 


Rout 


(ao) 


E v 


(au) 




Ev (cm" 1 ) 




T 

- 1 rev 


122 


92 


85,5 ps 


R-min 


-1.62 


x 10" 


-5 


-3.57 


196 


ps 


10 ns 


129 


98 


93.7 p 


-Rl 


- 1.21 


x 10" 


-5 


-2.65 


250 


ps 


15.3 ns 


137 


106 


107.48 p 


a Rmax 


- 0.79 


x 10" 


-5 


-1.74 


350 


ps 


15.7 ns 



TABLE III: Characteristic constants for three levels in the outer well of the 0^ (P3/2) potential 
curve, vtot is the vibrational number, while the numbering v is restricted to the levels in the outer 
well. — E v is the binding energy, R out is the outer turning point, T^b and T rev are respectively the 
classical vibrational period and the revival period as defined in Eqs. (J42I43|) in text. The level 



w=98 is identical to the (vq+2) level considered in Ref. 



■ El 



on tunneling. 



Due to the linear chirp parameter, an energy range of about 2H\x\tc = 1-74 cm" 1 is 
swept by the laser frequency, corresponding to 15 vibrational levels in the vicinity of v= 98, 
for which we have reported data. However, due to the spectral width of the pulse, levels in 
an energy range of 0.98 cm~ 1 on both sides can also be excited, and this will be analyzed 
with the numerical results. 



Since the coupling with the ground state may be involving the last least bound levels 
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in the ground state potential curve, we also report in Table IIVI the binding energies and 
vibrational periods for those levels, together with the time-delay for the continuum level at 
T =54.3 fiK, as defined in Eq. flSJ). 



V 


E (au) 


E (cm" 1 ) 


T V ib 


v"=50 


-3.87x 10~ 7 


£=-.084 cm" 1 


T vib =581 ps 


v"=51 


-1.23x 10~ 7 


£=-0.027 cm" 1 


r„ i6 =1460 ps 


v"=52 


-1.94x 


£=-0.042 cnr 1 


T vib =7.8 ns 


v"=53 


-2.43x 10" 11 


E=-5x 10 ~ 6 cm^ 1 




continuum 


+ 1.71 x 10" 10 


E=+3.77x 1Q- 5 cm- 1 


r ~ 476 ns 



TABLE IV: Constants for the last levels in the 3 £+(6s+6s) potentials. The present calculations 
are not including the hyperfine structure. 

The classical vibrational period has no meaning for the last level, since the wavefunction 
mainly extends in the classically forbidden region. We should note the very high value 
of the time delay in the present problem, due to the large value of the scattering length, 
demonstrates the strongly resonant character of the collision, and makes the collision time 
by far the largest characteristic time in the problem. Looking at Tables IIVI and [HJ we see 
that the spacing between the level w"=53 of the ground state and the continuum level with 
collision energy E = kT= 54.3 fiK « 3.77 x 10~ 5 cm _1 can be considered as negligible at the 
scale of the energy uncertainty imposed by ~ 30 ns radiative lifetime of the excited state. 



B. Description of the initial state 

The initial continuum state is represented by a stationary wavefunction in the ground 
state potential curve and describing the collision of two cold Cs atoms at the energy kT 
corresponding to T = 54.3 fiK . It is chosen with a node at the external boundary Lr of 
the spatial grid (see below). This wavefunction is drawn in Fig. EJ together with a typical 
gaussian wavepacket, showing how unlocalized is the initial state in the present work. This 



is an important modification compared to previous calculations 



|36( using a gaussian 



wavepacket: the population transfer can occur in a very wide range of internuclear distances, 
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well outside the resonance window, and in particular at large internuclear distances where 
a large density probability is localized. No velocity distribution is considered in the present 
work, this effect being treated in paper II. 



C. Numerical methods 



Details on the numerical calculations will also be given in paper II: they involve mapped 



grid methods 



3811 rather 



to represent the radial wavefunctions, using a sine expansion 
than the usual Fourier expansion, in order to avoid the occurrence of ghost states. The 
introduction of adaptive coordinates is necessary to implement a spatial grid with few points 
(1023), but with a large extension Lr = 19250 a® allowing representation of initial continuum 
states at ultra-low collision energy. Due to the large size of the grid, and the very small 
kinetic energy of the problem, it is not necessary to put an absorbing boundary condition at 
the edge of the grid. Since the initial state wavefunction is normalized in a box, the results 
given below will be depending on the value chosen for Lr. 

The time-dependent Schrodinger equation is solved by expanding the evolution operator 
exp[— i~H.t/h} in Chebyschev polynomia 4sj|. The time propagation is realized by discrete 
steps with a time increase much shorter than the characteristic times of the problem, already 
discussed in Sec. IIII Al and reported in the Tables ITTllI VI In the present problem, the smallest 
time scale is the local Rabi period at small distances, ~ 0.18 ps, which controls the time step 
At. In typical calculations, we have chosen At « 0.05 ps, the quantity H\l/ being evaluated 
112 times at each time step. 

The dynamics of propagation of the wavepackets in the ground and excited potentials is 
analyzed by studying the evolution of the population in both surfaces: 

P e (t) =< V e (R,t)\itf e (R,t) >, Pg(t) =< Vg(R,t)\V g (R,t) > (45) 

A more detailed information is provided by the decomposition of the wavepackets on the 
unperturbed vibrationals states v of both potential surfaces S = a 3 £+ or 0~: 

P Sv (t) = | < 9 Sv (R)\iS/ etg {R,t) > | 2 (46) 



18 



D. Results 



All the numerical results presented below are obtained from an initial state wavefunction 
normalized to unit in a large box of size L#=19250 clq. They correspond to branching ratios 
of the photoassociation process towards different final vibrational levels, eigenstates of the 
Hamiltonian for the unperturbed molecule. Due to the linear character of the Schrodinger 
equation, the population transferred to a given level is proportional to the normalization 
factor for the initial state. Therefore, in order to deduce the populations corresponding 
to an energy-normalized initial continuum state, for E— 1.71 x 10 -10 au (T=54 //K), the 
computed populations should be multiplied by the density of states dnj dE=l/ (8.87 X 1CT 12 ). 
Besides, since the calculations concern only one pair of atoms, and are not velocity averaged, 
the values for P e (t) should be scaled to get an estimation of the population transfer due to 
a single pulse in a trap containing N atoms: this will be done in Section I VI II 
The results of the calculations for the time-dependence of the populations are presented in 
Fig. both for a negative and a positive linear chirp parameter. The main conclusions are: 

• There is a population transfer from the ground 3 £+ state to the excited 0~ state which 
decreases by two orders of magnitude from P e (tp) ~0.032 near the maximum of the 
pulse, to P e (t >tp + to) = 3.2x 1CT 4 at the end of the pulse (see Fig. Ek,b). 

• Calculations considering the same laser pulse with an initial state represented by a 
gaussian wavepacket normalized to unity yield a much larger value (~ 0.59) for the 
final population in the excited state. However, the ratios between this final population 
and the probability density in the range 85 < R < 110 ao are similar (0.69 for the 
gaussian wavepacket, 0.73 for the continuum eigenstate). This result is a signature for 
the predominant contribution of the resonance window. 

• The maximum of the transferred population occurs ~ 5 ps before the maximum of 
the pulse for a negative chirp, and ~ 5 ps after this maximum for a positive chirp (see 
Fig. Ekb). 

• After the pulse, most of the population in the ground state is going back to the 
initial continuum state, but a small fraction 3 x 10~ 4 ) is transferred to the last 
vibrational levels in the 3 £+ potential, v"=5l-53 ( Fig. 01). This remaining population 
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is independent of the sign of the chirp. The population transfers towards bound levels 
of the excited 0~ and of the ground 3 E+ state are comparable. 

• While many bound levels in the excited 0~ potential are excited during the pulse, 
essentially the levels from v= 92 till v= 106 remain populated after the pulse (see Fig. 
Eh)- They correspond to the range of energy around v= 98 swept by the laser during 
the time window and presented above as a "resonance window" . 

• In contrast, highly excited levels from v=122 to v = 137 are significantly populated 
during the pulse (see Fig. Eh), but the population vanishes at the end of the pulse. 

• Whereas the oscillatory pattern during the time window markedly depends upon the 
sign of the chirp, the final population in the excited state, as well as the vibrational 
distribution, is nearly independent upon this sign. 

Such results demonstrate the existence of a "photoassociation windonF including all the 
levels in the energy range between v— 92 and v— 106 where population transfer is taking 
place. We shall provide interpretation for this result in the following section. 



V. A TWO-STATE MODEL FOR ADIABATIC POPULATION INVERSION. 
A. Impulsive approximation. 



For all the following developments, we shall use the impulsive approximation [33], as- 
suming that the relative motion of the two nuclei is frozen during the pulse duration. The 
kinetic energy operator being neglected in Eq. (j27jl . the two-level Hamiltonian becomes 

fi-H-f-ff-fr °) + ) (47) 

\0 V ) \W(t) -A(R,t) j 
where the first term introduces a /^-dependent phase while the dynamics is contained in the 
second term. 



B. The adiabatic basis 

When the impulsive approximation is valid, diagonalization of the hamiltonian H^i?, t) 
at each distance R will define a new representation, in the framework of a standard radiation- 
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driven two- level system |42j, |49j. In Eqs. (j47|) . the notations W and A are similar to the ones 
used in the textbook by Cohen- Tannoudji et al to describe a two-state model. Defining 
a time-dependent local Rabi frequency 



htt(R,t) = jA 2 (R,t) + W 2 (t), 



diagonalization of H 1 yields two eigenenergies E + and E_ 



E±(R,t) = V(R,t)±Ml(R,t) 



(48) 



(49) 



with two eigenf unctions, hereafter referred to as adiabatic, J 7 ^, J 7 ^. 

The set of two adiabatic functions is deduced from the diabatic functions ^ by a 

n 

rotation 1Z [50J: 



K 



\7L 



cos 2 sin 2 



— sm 2 cos 2 



the angle 9 being defined by the relations: 

\W(t)\ 



sin[6(R,t)\ 



™s[e(R,t)) = tm[e{R,t)\ - lw{t)l 



(50) 



htt{R,ty m(R,ty A(R,t) 

The two-channel wavefunction representing the evolution of the system can then be writ- 
ten in the adiabatic representation as 



\V{R,t) >= a + {R,t)\F% > +a~{R,t)\J^ > 



where a + (R,t) and a (R,t) verify 
/ 



a + (R,t) 
\a~(R,t) 



E+(R,t) 




a + (R,t) ihde{R,t) 1 
E„(R,t) I \ a-(R,t) ) 2 dt [-10 



(51) 



a + (R,t) 



\a-(R,t) 
(52) 

We shall discuss now adiabatic evolution when the second term in the r.h.s. of Eq.(|52|) is 
negligible. The validity of such an hypothesis will be discussed below in Sec JVI Al 

C. Condition for population inversion during the pulse duration. 

We assume that the effect of the pulse is negligible outside the time interval tp — 7> < 
t < tp + Tp, without assumption on the value of tf. Considering only the first term in the 
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r.h.s. of Eq. (|52|) . the solution of Eq. (J52j) becomes straighforward : assuming that before 
the beginning of the pulse, at t — tp — r F , there is no population in the excited state, the 
evolution of the diabatic wavefunctions is described by 



6(R,t) 
2 



B(R,t) 
2 



sin ») cos am 











™ J t p — Tp 



E- {R,t')d£ 



J 



X 



cos g&Sjpj sin 

-sin W-^) cos^f^ 







*»(R,tp-T F ) 



Introducing the accumulated angle 

a(R,t) = ~ [ E+(R,t')dt' 



'tp—Tp 

the two diabatic wavefunctions at time t are 



W2(R,t) = -[isin 



6t + 9i . . . . $t — 
sm a(it, rj + sin 



cosa(.R,t)] V"(R,t P -T F ) 



2 v ' ' 2 

Q \ Q. Q ^. 

\&"(R,t) = [i cos ' - sin a(fl, t) + cos * * cosa{R, t)} V"(R, t P - t f ), 
y 2 2 y 

where we have introduced the time- dependent angles 



(53) 



(54) 



(55) 



6 t , = 6(R,tp-T F ), 
t = 0(R,t). 



(56) 
(57) 



1. Pulse of finite duration 2 r F 

We are interested in the conditions leading to population transfer, at internuclear distance 
R, from the ground to the excited state, for a pulse starting at t = tp — r F and stopping at 
t = tp + r F : 

W{t P - t f ) = W{t P + 7>) « 0. (58) 
From Eqs.([5U|). we see that the rotation angles are such that 

sin 6(R, t P - r F ) = sin 6(R, t P + t f ) = (59) 

/)/ p , v A(R,t P -r F ) A(R,t P + r F ) 

cos6{R,t P -T F ) = -; cos9{R,t P + T F ) = — - 60 

\A(R,tp-T F )\ \A(R,tp + T F )\ 

According to the sign of the time- dependent energy-splitting, the angles 9(t F — t f ) and 
6{tp + r F ) can take the values or n. When the two angles have the same value, the 
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correspondance between diabatic and adiabatic states is the same before and after the pulse, 
i.e \^ e (R,t P + t f )\ = \W£ >e (R,t P - r F )\. In contrast, when the sign of A(R,t) defined in 
Eq. changes during the pulse, the two states are reversed. This population inversion 
takes place at a distance R provided that the level splitting between the dressed potentials 
verifies 

2\A L (R)\<h\ X \T F . (61) 

Then, 

\0{R,t P + T F )-0{R,tp-T F )\=n=> \^(R,tp + T F ) = \^ g (R,tp-r F )\. (62) 

For a given pulse (fixed values for x an d t f ), the condition (pTTj) for population inversion 
shows that the transfer will take place in a region of internuclear distances around Rl (such 
as Al(Rl) = 0), the extension of which depends upon the difference 2|Ax(i2)| between the 
two potentials and upon the energy range h\x\ T F swept during the pulse. 

In most photoassociation experiments the lasers are tuned so that Ri is located at large 
distances, where the difference 2\Al(R)\ between the ground and excited potential curves is 
very small, so that we may predict a large photoassociation window. 

2. Extension to a gaussian pulse. 

As it was shown in Sec III B 2l when the chirped pulse is gaussian, most of the intensity is 
carried during the time window (t — t F ) G [—t c , +tc], which suggests that the population 
inversion is mainly realized during this temporal window. Therefore, the condition for 
population inversion during the time window is written introducing tc instead of r F in the 
relation (JBT]) . 

• For internuclear distances R such that A(R, tp — r c ) > 0, the angle 0i(R,tp — tc) ~ 0, 
and the diabatic functions are 

WZ(R, t) = - sin ^e +ia{R ' t] ^(R, t P - t c ) 
2 y 

* u g (R, t) = cos |e +ia ^^( J R, t P - to) (63) 
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• In contrast, when A(R,tp — r c ) < 0, 9i(R,tp — tc) ~ 7r, and 



*" e (R, t) = cos -je~ ia ^^ g (R, t P - to) 
V^R, t) = sin 6 je- ia ^^ g (R, t P - r c ) 



(64) 



Let us emphasize that in the adiabatic following of an instantaneous eigenstate of the 
system, since the angle 0j stays strictly equal either to or ir, the population in the ground or 
in the excited state does not depend upon the accumulated phase a(R,t). This is not valid 
when sin^j ^ 0, cos#j ^ ±1, since the population at a given R value in the excited channel 
exhibits Rabi oscillations with a time period corresponding to a tt variation of the accu- 
mulated angle a(R, t). The latter situation requires that the laser field is turned on suddenly. 

The population inversion condition defines a range of distances where photoassoci- 
ation is taking place, equivalent to the resonance window defined above in Eqs. 
However, the validity of such an interpretation relies upon the validity of the adiabatic ap- 
proximation both inside and outside the window, and of the choice for the temporal window. 

The conclusions are independent of the sign of the chirp, which is no longer true when 
non-adiabatic effects are involved. 

VI. ADIABATICITY REGIME IN THE PHOTOASSOCIATION WINDOW AND 
CONDITION TO AVOID RABI CYCLING AT LONG RANGE 

A. The adiabaticity condition 

From Eq. (j52j) . we see that the system will present an adiabatic evolution provided the 
non-diagonal term is negligible compared to the energy difference (E + — E_): 



From Eq. (JHOJ) we may derive 89/ dt = [A(dW/dt) - W(dA/ dt)}/[(hVL) 2 }). From Eqs. 
f)28|29jl . the time-derivatives of the coupling term W(t) and of the level-splitting A(R, t) are 
easily derived, leading to the explicit form of the condition (J65J) . 



89 
dt 



<4fi( J R,t) 



(65) 



n 



^(t-t P )A(R,t) + 




h X . _4[(A( J R,t)) 2 + (H/(t)) 2 ] 3 / 2 



4 



[hn(R,t)] 3 



(66) 



c 
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that we wrote as an inequality between two quantities with the dimension of an energy- 
square. The quantity in the left-hand side becomes large in the wings of the pulse (when 
dW/dt is large), far from resonance (when A(R, t)) is large) and for a large value of the chirp 
rate x- The quantity in the right-hand side is the smallest at the instantaneous crossing 
point R c {t) such that A(R c ,t) = 0, where htt(R c ,t) = W{t). 
The condition (|66|) takes simple forms if one chooses particular values of R or t: 



1. At t = tp , when the coupling reaches the maximum W max = Wl^j^, one can write 
an adiabaticity condition function of R: 

fi \x\ < , — — — (67) 

W l Jt l /t c 



2. In particular, when both t = tp and A L (R L ) = 0, the condition at the "crossing point" 
Rl for the two dressed potentials, when the coupling reaches its maximum value, reads 

h\ X \r c « 8W 2 L ^. (68) 

3. It is also worhtwhile to consider the very strict adiabaticity condition at the 
instantaneous crossing point Rc{t)'- 

h 2 \ x \ <8(iy(t)) 2 . (69) 

The condition (}6T?|) will not be verified in the wings of a gaussian pulse. Nevertheless, 
if we consider the time window \t — tp\ < tq already defined above, the coupling 
parameter has the lower bound \W max , so that during the time window the condition 
is simply: 

Mxlrc « Wl^-. (70) 
The validity of such hypothesis will be discussed below. 



B. Validity range of the adiabaticity condition 

We shall first ensure that, during the temporal window, the adiabaticity condition is 
indeed valid all across the photoassociation window, defined both by Eq. (J31|) and by the 
relation (|61|) in which Tp = to- Then we shall check that no population transfer is occuring 
outside the window, therefore justifying the definition of a window. 
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1. Adiabaticity condition within the photoassociation window 



From Eq.fpj)). we define an adiabaticity parameter, which should be -C 1, as 
Y(BA-Z»- ^/rc)\-^2(t-t P )/(r c )A(R,t)+hxrc/2\ 

1 ' j x d 4[hn(R,t)] 3 /w(t) • [n) 

During the time window, and within the photoassociation window, the local instantaneous 
level splitting verifies \A(R,t)\ < \hxrc\, so that the numerator X n has an upper limit 
X n < fr 2 |x|(41n2 + 1/2), reached at the edges of both the time window and the distance 
window. Under the same conditions, a lower limit of the denominator can be found from 
the pulse shape, considering both the upper and the lower limits on the coupling during the 
time window, defined above in Eq. (|24|). Therefore, a sufficient condition for adiabaticity 
within the photoassociation window is 



I ^(4hi2 + l/2) 
16 x {h\x\Tc) = 52.36ft \x 



1 + ,fj C ^ < w l => X(R,t) < 1. (72) 



t l " ' IA '^ (4 In 2 

The condition (fT2|) means that the laser intensity is sufficient to yield a coupling much larger 
than the geometric average between the energy range swept by the pulse {2Ti\x\tc) and the 
energy width of the power spectrum |£(u;)| 2 , which is (Mln2/ri). The same condition may 
be written as an upper limit on the energy swept by the chirped pulse through 

52.36ft 2 |x|r c < W 2 L r L} (73) 

where in the rhs of Eq. (|T3*j) we recognize the quantity WItl, proportional to the energy 
carried by the field. This energy should be sufficient to cause adiabatic transfer. 
Finally, it is instructive to rewrite the condition (|73*jl in terms of the time scales T^ abi and 
T c hir P , introduced above in Eqs. (j34J37|) . as well as the width tl or the time scale T spect 
associated to the spectral width of the pulse through Eq. (j3*%|) : 



TroM 0. 24 ^T c hir p TL — 0.16 \jT c hirpT spec t. (74) 

Optimizing the adiabaticity condition requires the increase of the intensity of the pulse or 
the reduction either of the photoassociation window (large value of T c h irp ) or of the spectral 
width of the pulse (large value of tl or T spect ). Since the population transfer will depend 
upon the width of the photoassociation window, a compromise has to be chosen. 
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2. Adiabaticity in the asymptotic region 



It is often convenient to ensure that no population is transferred outside the photoasso- 
ciation window, in the region of distances where the two dressed curves never cross during 
the time interval (t — t P ) G [—r c , +tq}- In particular, we have considered the long distance 
region, where 2A^(i? — > oo) = 5£* being the detuning relative to the atomic transi- 
tion, as it is defined in Eq. (fT|). In this region, considering upper and lower limits of the 
coupling and of the level splitting, as in the previous section, the validity of the adiabatic 
approximation is ensured under the condition 

= 5?5=[(? - + (^p)r 3 '^25f + " l)t|Tc(4 '" 2+1) | « 1. (75) 

4tc 2 2 4 " 2 

This condition can be simplified under some circumstances, which depend upon the rel- 
ative values of W max , and ^|x|rc. 

VII. DISCUSSION 

From the simple model developped above, we shall analyse further the numerical results 
of Sec. HV1 

A. Analysis of the numerical results in the framework of the two-state model 

First, we have represented in Fig|H] the variation of the local frequency Q(R,t), defined 
in Eq . (J48|) . as a function of R for three choices of the time t = tp,tp ± Tc- The long 
range minimum in Q(R,t) moves from R max = 110 ao to R m in =85 ao during the temporal 
window [—35, +35] ps. Such minima occur at the distances where the local time-dependent 
energy splitting A(R,t) goes through 0, as illustrated in FigEl and therefore define the 
photoassociation window. Further minima in the local frequency occur at shorter range, 
corresponding to inner crossings of the two dressed curves which do not play a role in 
the photoassociation process, since the amplitude of the initial wavefunction \ty init (R)\ 
is negligible at such distances. Also indicated in FigEk are the values of the local Rabi 
period at time t = tp, further illustrating the wide range of variation with R, from 12 ps 
at infinity to 34 ps at R L , with very small values (0.18 ps) and larger ones (28 ps) at short 
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range where the two curves part and cross again. This behaviour is linked to the particular 
shape of the potential curves, with a double well behaviour for the excited one. Looking at 
Fig. Eh, it is clear that due to the reduction of the pulse intensity in the wings, reducing 
W(t) by a factor of 4 at tp ± tc, the minima in Q(R,tp ± tc) are deeper at the borders 
of the time window, multiplying the Rabi frequency by a factor of 4. This is why it will 
be difficult to verify the adiabaticity condition at the borders of the time window and the 
photoassociation window. 

When the adiabatic approximation is valid, the impulsive two-state model described in 
Sec. predicts population inversion within a photoassociation window characterized by a 
change of n between the angle 9(R,tp — t c ) and the angle 9(R,tp + t c ), where 9(R,t) is 
defined in Eq. ()5(Jj) . The variation of the angle 9(R,t) as a function of R is represented in 
Figdfor the same three choices of the time. We see that in a range of distances extending 
from 85 to 110 ao, 9 varies approximately by it from tp — Tc till tp + r^: a population 
inversion can indeed be predicted. 

The rather good agreement between numerical calculations and the predictions of a simple 
two-state model seems to indicate that the adiabatic approximation is indeed valid in our 
problem. In order to get more insight, we have studied the variations of the adiabaticity 
parameter X(R,t) defined above in Eq l71l which are represented in Fig. |HJ 
The adiabaticity parameter takes maximum values at the border of the time window (t = 
tp ± Tc) and of the photoassociation window R = R m in, Rmax, as illustrated in FiglHl^a). 
For the pulse with negative chirp described here, we find two adiabatic regions and one 
non-adiabatic: 

• The adiabaticity condition is verified within the photoassociation window, as illus- 
trated in Fig. IHl^b); 

• It is only approximately verified at R ~ R m in for t ~ (tp + tq) and at R ~ R ma x for 
t ~ (tp — tq), where \X(R=100 ao,tp — tc)\ reaches the value 0.5, as illustrated in 
FigHb); 

• It is not verified out of the photoassociation window, in the vicinity of R ma x and R m in, 
at distances i?=110 and 85 ao as illustrated in FigJH^c); 
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• In contrast, the adiabaticity condition is well verified in the asymptotic region, as illus- 
trated in FiglHfd) . We should note however that the present pulse has been optimized 
for that purpose. Other choices yield an important population population transfer at 
large distances. 

The present discussion of the adiabaticity criteria can be attenuated when using a nar- 
rower time window, our choice of [tp — tc, tp + tc], corresponding to transfer of as much as 
98 % of the energy in the pulse, being probably too severe. A sufficient choice considers the 
time window \t — tp\ < O.Qtc, limited by the two inflexion points in the coupling term W(t), 
during which 84% of the intensity is transferred. Better boundaries for the adiabaticity 
parameter are then obtained: 

\t - t P \ < t c -> 1.99 x 1(T 3 < \X\ < 29.99 
\t - t P \ < 0.6t c -> 13.3 x 1(T 3 < \X\ < 1.3849 (76) 

B. Discussion 

It is instructive at this stage to refine our interpretation of the numerical results 
presented above in FigJSJ Due to the coupling by the laser, vibrational levels in the 0~(P3/2) 
are excited. We shall call "resonant levels" the t>=92-106 levels, corresponding to the 
range of energy swept by the instantaneous frequency u(t) during the time window, and 
"off-resonance levels" the other ones, corresponding to energies swept by the wings of the 
spectral distribution during the time window, or by the whole spectral distribution outside 
of the time window. Among them, excited levels are numerous, such as the t>=122-137 
range, due to the high density of levels close to the dissociation limit, and to the good 
Franck-Condon overlap with the initial state. 

When the chirp is negative, the instantaneous crossing point is moving from R ma x to R m in, 
so that the most excited levels are populated first. This is illustrated for the "resonant 
levels" in FigjSfc). In contrast, for a positive chirp, since the resonance window is swept 
from R m in to R m ax, the lower levels are populated first, as illustrated for the resonant levels 
in FigEtf). 

For the "off resonance" levels in the excited state, a maximum in the population P v {t) 
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(see Eq. |46J) appears, independently of the level, 4.45 ps before tp for x < 0, and 4.45 ps 
after tp for x > (FigJ^e)). A quantitative evaluation of this advance (resp. delay) can 
be obtained from analysis of the adiabatic population transfer at large distance R > 150 
a , where 2A L (R) « the adiabatic parameter being such that \X(R,t)\ < 0.03 (FigJHJ). 
The diabatic wavefunction in the excited potential 0~ is given by Eq.(|63|) where both the 
angle 8(R,t) and the accumulated angle a(R,t) are independent of R. Furthermore the 
vibrational wavefunction in the excited potential curve is strongly localized in the asymptotic 
region. Therefore the time-dependence of the population P e (t) as defined in Eq.pS^I scales as 
sin 2 |(oo,t), reaching a maximum for a time t m satisfying ^|(oo,t m )=0, and corresponding 
to a zero value of the adiabaticity parameter. Trivial calculations yield the condition 



tm ~ tp ~ 4^% (77) 



which gives -4.37 ps for the present pulse with negative chirp, and explains the symmetry 
when changing to a positive chirp. The excellent agreement between this estimation and 
the ±4.45 ps value found in numerical calculations is a signature of the validity both of 
the adiabatic approximation and of the impulsive approximation in the asymptotic region. 
Indeed, the vibrational period for levels v > 122 is T V ib(v) > 780 ps. The maximum is also 
visible in the variation of total population in Fig|S^b), therefore dominated by asymptotic 
excitation. 



For the "resonant" levels, oscillations in P v (t) appear only for the levels which come 
to resonance with the instantaneous frequency uj(t) before the maximum of the pulse, at 
a time tc(v) such that tc(v) < tp. For x < 0, oscillations are therefore observed for the 
highest levels, like v—107 and 102 in FigJS^c). For x > 0, the population of the lowest 
levels, like v=95 in FigGJf) is oscillating. For the v=98 level, resonant at the maximum 
of the pulse, the time-dependence of the population does not depend upon the sign of the 
chirp (red curve in Fig|3Jc,f )), and oscillates weakly. Such oscillations are a signature of 
significant non-adiabatic effects in the population transfer at the instantaneous crossing 
point around t = tp — re- We have shown in Sec iVIl how the adiabaticity criterion is 
most difficult to verify at the instantaneous crossing point and in the wings of the pulse. 
Away from the adiabatic regime, the passage through resonance does not lead to a total 
population transfer: coherent excitation of a two-level system with linearly chirped pulses 
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results in "coherent transients" previously studied in the perturbative limit 
Such transients are governed by interferences between the population amplitude transferred 
at resonance t = tc(v), and population amplitude transferred after, for t > t c {y). The latter 
amplitude is significant provided the maximum of the pulse occurs after the resonance, 
so that tp > tc(v), resulting into the observed oscillations. The values of the observed 
periods can correspond to a Rabi period at the instantaneous crossing point: for v=98, 
resonant at t — tp, the period of 36 ps is close to TRab(RL,t p )— 34 ps; for the level v=102 
with a negative chirp, or v =95 with a positive chirp, the 25 ps period is associated to a it 
variation of the accumulated angle a(R,t) calculated at R distances corresponding to the 
outer turning points of these vibrational levels. 



Looking back to the derivation of the equations describing adiabatic transfer in SecfVl 
we see that in the case when sin 0j ^ ( or ir), i.e. when the adiabatic evolution starts only 
once the laser has been turned on for a while, the evolution of the population in the excited 
state is no longer described by Eq. (or Eq.(JBH)), but verifies 

9 9- 9 9 1 

\^(R, t)\ 2 = [sin 2 ^ cos 2 i + cos 2 ^ sin 2 ^ - - sin 9 t sin 0* cos[2«(i?, t)]\V"{R, t P - r ad \ 2 

(78) 

for x < 0, where we have called tp — r a d > tp — tq the time where the adiabatic evolution 
starts (a similar equation can be written for \ > 0, with change of sin into cos). 



The total population at a given R value in the excited channel therefore exhibits 
Rabi oscillations, with a time period corresponding to a 7r variation of the accumulated 
angle a(R,t). The modulation rates appearing in Figl^c) for the high resonant vibrational 
levels with large vibrational periods T v n, ~ 350 ps, which are the first populated with 
a negative chirp, are much larger than those observed in Fig|Sff) for the low resonant 
vibrational levels, with T v n, ~ 200 ps, the first populated with a positive chirp. Rabi 
oscillations being manifested in cases where the vibrational period is much larger than the 
Rabi period, a frailty in the impulsive approximation can be supposed for levels for which 
these two characteristic times become comparable. Besides, as discussed by Banin et al 
(1994) (33(1, the wavepacket created in the excited state being accelerated towards short 
internuclear distances, the breakdown of the impulsive approximation is more severe for 
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X < than for x > 0. Indeed, for x < 0, the instantaneous frequency of the laser remains 
resonant during the motion of the wavepacket and then recycles the population back to 
ground state |35| . 

Therefore, the simple model for two-state adiabatic population inversion, developped in the 
impulsive limit, allows a qualitative interpretation of the numerical results. In particular, 
the vibrational levels in the excited state where population remains transferred after the 
pulse are well predicted by this model. Also is the weak dependence of the final population 
as a function of the sign of the chirp. However, the phase of the probability amplitude on 
each level is not independent of the sign of the chirp, so that after the pulse the vibrational 
wavepacket evolves in a very different way for positive and for negative chirp. Furthermore, 
detailed analysis of the numerical results indicates the limits of both the impulsive model 
and the adiabatic population transfer. 

VIII. APPLICATION: CONTROLLING THE FORMATION OF ULTRACOLD 
MOLECULES VIA FOCUSSING AND IMPROVING PRESENT EXPERIMENTAL 
SCHEMES 

A. Fo cussed vibrational wavepacket 

The pulse described in Table HJ with negative linear chirp parameter x, has been designed 
in order to achieve focussing of the vibrational wavepacket at a time tp + where 
T V ib/2= 125 ps is half the vibrational period of the level v=98. The levels in the range 
from approximately t>=106 till v =92 are populated in turn, the chirp parameter necessary 
to compensate the dispersion in the vibrational period of the wavepacket being chosen as 
X = —2iiT rev /(T V ib) 3 , i.e. adjusted to match the revival period T rev (see Eq. (H3*|) ) of the 
resonant level t>=98. Indeed, for x satisfying Xv — ~ 2^T rev (vo)/[T vib (vQ)] 3 , the amplitude 
resonantly transferred at t = tc(v) creates a wavepacket at the outer turning point of each 
vibrational level v ~ Vq, which reaches the inner turning point at time t = tp + T vi b(vo) 
independently of v. The evolution of the wavepacket after the pulse has been computed, 
and we have represented in FigJH] a snapshot of the wavepacket in the excited state at time 
tp + T v n,/2. 

• An important focussing effect is visible for a negative chirp, the wavepacket presenting 
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a peak in the vicinity of the inner turning point of the v=98 level. 

This effect is reduced when a positive chirp is used, a large part of the population 
being in the large distance region, due to the late excitation of the upper levels v= 
98-107, with a much larger vibrational period. 

When there is no chirp, the population transfer is much smaller, the factor w 3 on the 
maximum amplitude resulting into one order of magnitude in the transfer probability 
( 4 xlO ~ 5 instead of 3.5 x 10 _4 ). Besides the wavepacket is no longer focussed since 
the levels that have been populated are no longer in the vicinity of v=98, but belong 
to the domain v > 103 of very excited levels; from the behaviour of this wavepacket 
it is clear, as illustrated in FigEfb), that the levels have a large vibrational period. 
The distribution of vibrational levels is governed by the spectral width hSu and by 
the overlap integral between the initial and final wavefunctions, which is dominant for 
v « 130. 

Focussing of the wavepacket on the inner classical turning point of the vibrational 
levels will allow to consider population transfer towards bound levels of the ground 
triplet state potential curve via a second laser, in a two colour pump-probe experiment. 
It should then be possible to populate efficiently low v" levels in the ground state. 

For such purpose, the optimization of the pulse can be achieved through analytical 
formula, considering the revival period, and the scaling laws governing the dynamics 



of long range molecules \\. 153]. 



B. Towards new experiments 

A proper estimation of the photoassociation rate requires incoherent average over a ther- 
mal distribution of energy normalized initial states. These energy-normalized continuum 
wavefunctions are deduced from wavefunctions normalized to unit in the box of size Lr, 
accounting for the density of states in this box. Furthermore, following Ref. 2f|, for low 
temperature, it is possible to estimate the photoassociation yield by assuming a weak vari- 
ation of the population as a function of the initial energy. For an assembly of atoms in a 
volume V= 5x 10~ 4 cm 3 , with density N at =10 n cm" 3 , the number of molecules formed in 
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one pulse under the conditions discussed in the present paper is 0.690. With a repetition 
rate of 10- Hz, the yield is 6.0x 10' molecules per second, which is significant iarger thai 
rates obtained with continuous lasers 141. 



IX. CONCLUSION 

The present theoretical paper has investigated the possibility offered by replacing con- 
tinuous lasers by chirped laser pulses in ultracold photoassociation experiments. We have 
performed and analyzed numerical calculations of the population transfer due to a chirped 
pulse, from a continuum state (T « 54fiK) of the ground triplet state Cs2 a 3 £+(6s + 6s), 
to high excited vibrational levels in the external well of the 0~ (6s + 6P3/2) potential. The 
central frequency of the pulse was chosen resonant with the level v —98 of the 0~ external 
well, bound by 2.65 cm -1 . We have used a pulse with gaussian envelope, having a peak 
intensity II ~ 120 kW cm -2 , a linear chirp parameter of -4.79 x 10~ 3 ps -2 , and temporal 
width 7~c = 34.8 ps. Due to the linear chirp parameter, a range of energy of 2 x 0.87 cm -1 
is swept by the laser frequency, corresponding to resonance with 15 vibrational levels in the 
vicinity of v= 98, and referred to as "resonance window". The spectral width of the pulse, 
which remains independent of the chirp, covers a typical range of energies of ~ 0.98 cm -1 on 
both sides. 

The initial state of two colliding ultracold cesium atoms is represented by a continuum 
stationary wavefunction, with a realistic scattering length and nodal structure. This is 
oossible owing to use of a mapped sine grid method recently developed by Willner et al 
3sj|. The time-dependent Schrodinger equation is solved through expansion in Chebyschev 
polynomia 

Due to the large value of the scattering length, the initial wavefunction has a large 
density probability at large distances, resulting into a large Franck-Condon overlap with 
vibrational levels in the excited state close to the dissociation limit. Numerical calculations 
show that, during the pulse, a large amount of population is indeed transferred to levels 
close to the CS2 (6s + 6P3/2) dissociation limit: however, with the particular choice of the 
chirp, we have shown that this population is going back to the ground state after the pulse. 
In contrast, the 15 levels within the photoassociation window are coherently populated, 
and part of the population remains after the pulse. The same conclusion is obtained by 
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changing the sign of the chirp. Besides, an interesting result is that, due to the coupling at 
large distances between the two electronic states, a strong population transfer to the last 
vibrational levels of the ground a 3 S+(6s + 6s) is realized, making stable molecules with a 
rate as important as for the photoassociation into the excited state. This result should be 
further explored. 

The interpretation of the strong population transfer towards the 15 levels in the resonance 
window, independent of the sign of the chirp, is given within the impulsive limit, in the 
framework of a two-state adiabatic population inversion model. The numerical results can 
be qualitatively and even quantitatively interpreted by this model, defining the conditions 
for full population transfer within a "photoassociation window" covering a range of distances 
85 ao < R < 110 ao, and no transfer outside. Therefore the chosen chirped pulse is capable 
of controlling the population transfer towards the excited state, ensuring that no transfer 
occurs outside a chosen energy range. We have discussed the general conditions for adiabatic 
population inversion during a "time window" [tp — Tc,tp + Tc] centered at the time tp of 
the maximum of the gaussian pulse, and corresponding to 98 % of the energy transferred by 
the pulse. The adiabaticity condition within the "photoassociation window" and "during 
the time window" requires that the coupling should be larger than a quantity proportional 
to the geometric average between the energy range swept by the instantaneous frequency of 
the chirped pulse and its spectral energy width. To keep adiabaticity, a compromise between 
increasing the laser intensity and decreasing the width of the photoassociation window must 
therefore be found, and it can be discussed within the model we propose. Even though in 
the present calculations the adiabaticity condition is verified within the photoassociation 
window, a few non adiabatic effects were identified and interpreted. Further work will 
estimate more thoroughly the lower limit for intensity. 

As an example of control, the present pulse has been chosen in order to achieve focussing 
of the excited vibrational wavepacket at the inner turning point of the v =98 level, at a time 
delay after the pulse maximum equal to half the classical vibrational period T v n,{v = 98). 
The relevant chirp parameter was optimised from the revival time. Such calculations suggest 
a scheme for a two colour experiment where a second pulse, delayed by would transfer 

photoassociation population to low v" levels of the ground triplet state, now forming stable 
molecules in low vibrational levels. 
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In a forthcoming paper, we shall give more quantitative results for comparison with 
experiment by averaging on the initial state velocity distribution to compute the photoasso- 
ciation rate for a series of pulses with a variable repetition rate, considering different choices 
of the final vibrational levels by varying the detuning, and different laser intensities. We 
shall propose various schemes to optimize the efficiency of the population transfer, taking 
advantage of the conclusion of the present paper where we have shown that it is possible 
to design chirped laser pulses so that the photoassociation reaction is restricted to a well 
defined window. 

From the present calculations it is already possible to conclude that photoassociation 
with a chirped laser pulse is a very promising scheme. In the present example where the 
peak intensity of the laser was chosen II ~ 120 kW cm -2 , the estimated rate of ~ 7 x 10 7 
molecules per second in typical experimental conditions is definitely larger than what has 
been computed and measured for cw lasers. Since the present calculations were performed 
in a frequency domain close to a photoassociation minimum, much higher rates are to be 
expected by varying the detuning, and this will be the subject of a forthcoming paper. 
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FIG. 1: Scheme of the photoassociation process with a negative chirped pulse considered in the 
present work, illustrated in the case of Cs2- The potentials curves correspond to the ground triplet 
state a 3 S+(6s + 6s) and to 0~ (6s + 6p 3 / 2 ) excited electronic state. In the present work, the 
energy of the initial continuum state is neglected in the definition of the resonance condition. The 
double well behaviour in the excited curve is a particular feature of the chosen symmetry. 
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FIG. 2: Properties of a chirped pulse, (a) Stretching of the temporal width from tl to tq '■ 

variation of the amplitude for a transform- limited pulse of width tl = 15 ps ( solid line), and for 

the chirped pulse (broken line) obtained with linear chirp \\\ = 170 ps -2 . Note the reduction of the 

intensity at maximum, as described in Ea.l|21|) in text, and the broadening to tq = 34.8 ps. The 

time window [tp — Tc,tp + tc] is indicated by the horizontal line, (b) Central frequency sweeping: 

assuming that the carrier frequency ujl at t = tp is at resonance with the vibrational level v = 98 

in the Cs2 0~ external well potential, we have indicated the variation of ui(t) (see Ed. 11188 in text) 
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around U£ by %{p ~ tp)- 
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FIG. 3: Wavefunction of the system before the pulse; the range of distances [R m in = 85 ao, R ma x = 
110 ao], later on referred to as "photoassociation window" is indicated on the upper horizontal scale, 
the crossing point i?L=93.7 ao on the lower one. 

a) Variation of 10 x \Wjj L (R, 0)| (broken line), the continuum wavefunction describing the relative 
motion of two Cs atoms in the potential a 3 S^ at energy E = kT = 1.72 x 10~ 10 au, T = 54/iK. 
For visibility, this wavefunction, normalized to unit in a box of size 19 250 ao, is multiplied by 
10. An energy-normalized function can be deduced by dividing the unit normalized function by 
^JdE/dn = V8.87 x 10~ 12 . Also represented is a gaussian wavepacket of width a = 15 au, centered 
at R g = 95 ao ~ Rl- 

b) and c) Integrated density of probability J^\^\ 2 dR for the localized gaussian wavepacket (b) 
and the delocalized continuum wavefunction (c) : note that in the region 85 ao < R < 110 ao, the 




FIG. 4: a) Variation of the local time-dependent half energy splitting A(i?, t), defined in Eq, (|29|) as 
a function of R, for the time values t = tp ( solid line), at the maximum of the pulse; t = tp — tq, 
(dash-dotted line) at the beginning of the time window; t = tp + rc, (broken line) at the end of the 
time window. The distances where A(R, t) = define the crossing points Rl, Rcitp + t~c) ~ Rmin 
and Rc(tp — tq) ~ Rmax- b) Photoassociation window: in the reflexion model, the levels v=92 to 
v=106 can be populated when the crossing point Rc(t) is sweeping the range [R m in, Rmax]- The 
two vertical lines at short R correspond to a local maximum, related to the particular double well 
structure in the excited curve, having a potential barrier. 



44 



0.9B5 



0,975 



0.965 



" total "~^\ 


ft — - 


■ \ I 






/ / xs0 - 


1 . 





T 



i 1 1 1 

f*Wt 



black Z<0 /' r « 
blue *>0 // 



i.Utf-LK 



.5e-D4 



1 0e-04 



5.0e-D5 



0.0e+00 



50 100 150 200 250 50 100 150 200 250 

4.0e-D4 



030 



020 



010 



0.000 



I 'It'-.' 1 




total " / / 

: 1 








■(b) Jpr 1 





v=t)5-L07 




3.Qe-D4 



2.0e-D4 



- 1;0e-D4 



50 100 



1 50 200 



250 50 



O.De+DO 



e.Oe^OS 



4,0e-05 



2 0e-D5 



O.Oe+00 




100 150 200 250 

e.0e-05 



4 0e-05 



50 100 



150 200 250 50 

tips! 




150 200 

t(ps) 



2.0e-05 



D.Qe±OQ 



FIG. 5: Time- variation of the computed populations for the ground 3 S+ and excited 0~ electronic states 

during the pulse duration. The pulse is centered at tp — 150 ps. Black lines in Fig. (a), (b), (d), (e) 

correspond to negative chirp, and blue lines to positive chirp, (a) Full lines: variation of the total population 

in the ground state. Note that the minimum occurs w 5 ps earlier than tp = 150 ps with negative 

chirp, and « 5 ps later with positive chirp. Broken lines: variation of the population in the initial continuum 

level, at £7=1.72 x 10~ 10 au. (b) Variation of the total population in the excited 0~ state for negative (dark) 

and positive chirp (blue), (d) Population of the three last bound levels of the 3 E+ state: full line v"=51, 

long-dashed line u"=52, dash-dotted line u"=53. (e) Variation of the partial sums of populations on the 

bound vibrational levels in the 0~ potential: full line for low excited levels, v=95 to 107, resonanly excited 

during the temporal window; dashed line: highly excited levels, v— 122 to 135, off-resonance. Black and 

blue correspond to negative and positive chirp, (c) Details on the variation of the individual populations 
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of the g levels resonantly excited during the temporal window, for negative chirp: in particular red curve, 





FIG. 6: Variation of the local Rabi frequency Q(R,t), denned in Eg . (|48|) and corresponding to the 
pulse defined in Table [H as a function of the distance R, for three values of time, (a) Solid line: 
at t = tp, at the maximum of the pulse. Also indicated are the values of the corresponding local 
Rabi period Tp a b, see Ea. (|36j) in text, (b) Dash-dotted line: t = tp — tq, at the beginning of the 
time window. Broken line : t = tp + tq, at the end of the time window. Also indicated are the 
values of the characteristic times T v n, and T^ abi (tp,tp ± tc) significant for the vibrational levels 
resonantly excited at the moments tp + Tc,tp,tp — re, which have the outer turning points in 
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FIG. 7: Population inversion and photoassociation window, a) Variation of the angle 8(R,t), 
denned in Ea. (|5U|) for the same three time values as in Fig. HO The domain where 8(R, t — tq) and 
0(R, t + rc) differ by n is indicated by an horizontal arrow: it corresponds to population inversion, 
as discussed in Section IV CI 
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FIG. 8: Interpretation of the calculations: validity of the adiabatic approximation within the pho- 
toassociation window. The adiabaticity parameter X(R,t), defined in Ea.(|71jl in text, is represented 
for various values of time and distance. 

(a) Variation of X(R, t) as a function of R for t = tp (solid line), t = tp — tc (dash-dotted curve), 
t = tp + tq (broken curve) . 

(b) Adiabatic behaviour within the photoassociation window : Variation of X(R, t) as a function 
of time for R = Rl (solid curve), R = 90 ao (broken curve), R = 100 ao (dash-dotted curve). Note 
that maximum values are reached for \t — tp\ = tc- 

(c) Non-adiabatic behaviour at the border of the window: same as FigJHJa for R=85 ao (broken 
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curve) and R=1W ao (dash-dotted curve). 
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FIG. 9: Optimization of the chirp in view of focussing. Variation of the amplitude of the wavefunc- 
tion in the excited state \^f exc (R, t = tp + 125) | , at a delay = 125 ps after the maximum of 

the pulse (i.e. well out of the time window), where T v ^ is the classical vibrational period for the 
resonant level t>=98. (a) negative chirp, for the pulse defined in Tabled (b) same with x = 0j ( c ) 
same with positive chirp parameter. 



